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A set of relativistic mean field models is constructed including the Hartree and Hartree-Fock 
approximation accounting for the exchange of isoscalar and isovector mesons as well as the pion. 
Density dependent coupling functions are determined to reproduce the components of the nucleon 
self-energy at the Fermi surface, obtained within the Dirac-Brueckner-Hartree-Fock (DBHF) ap- 
proach using a realistic nucleon-nucleon interaction. It is investigated, to which extend the various 
mean field models can reproduce the DBHF results for the momentum dependence of the self- 
energies and the total energy of infinite matter. The mean field models are also used to evaluate 
the bulk properties of spherical closed-shell nuclei. We find that the Hartree-Fock model allowing 
for the exchange of a, u, p, 8 mesons and pions, yield the best reproduction of the DBHF results 
in infinite matter and also provides a good description of the properties of finite nuclei without any 
' adjustment of parameters. 
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I. INTRODUCTION 



One of the main challenges of theoretical nuclear physics is the attempt to develop a microscopic theory, i.e. without 
any adjustment of parameters, for the description of bulk properties of finite nuclei, which is based on high precision 
free space nucleon-nucleon (NN) interactions. This means that this approach connects properties of baryons in the 
vacuum (NN scattering data) with nuclear systems at densities around the saturation density of nuclear matter. 
Therefore, this approach should have a high predictive power, when it is used in extreme cases such as in highly 
isospin asymmetric nuclear systems. The study of these exotic nuclear systems is a fast-growing field of physics. 
The great scientific potential is demonstrated by the plans for large-scale exotic-nuclear-beam facilities, such as the 
the future GSI facility in Germany, and is driven by the expectation of observing nuclear properties, which are very 
different from those encountered so far, i.e. near the valley of stability. 

However, most of these attempts have already failed in completing the first milestone, which is the microscopic 
^{■^ description of the saturation properties of infinite nuclear matter in terms of such free space NN interactions^- 4 1 . 
Therefore three- nucleon forces have been introduced to fit the saturation point and/or properties of light nuclei [5|, |6|. 
Relativistic many-body approaches, in particular the Dirac-Brueckner-Hartree-Fock (DBHF) approach, however, have 
been successful in describing the saturation properties of nuclear matter Q-QjJ without the necessity to introduce 
many-body forces. Therefore a significant part of the 3-body terms, which are required in non-relativistic investi- 
gations, may represent the relativistic effects originating from the so-called Z-graphs in the expansion of relativistic 
propagators [3 HU . 

Although DBHF calculations have been quite successful in describing nuclear matter, full Dirac Brueckner cal- 
culations are still too complex to allow an application to finite nuclei at present, since the consistent treatment of 
correlation and relativistic effects for finite systems is a rather involved problem. In fact, Van Giai et al. [2(| ad- 
dressed this problem as one of the main open problems in Nuclear Physics. Different approximation schemes have 
been developed, which treat either the relativistic effects or the correlation effects in an approximative way. 

In the former approximation scheme, the Dirac effects are treated via a kind of local density approximation (LDA), 
whereas the correlations effects are taken into account by solving the Bethe-Goldstone equation directly for the finite 
nucleus under consideration l21H23j . This means that the self-consistency requirements of a conventional BHF 
calculation for a finite nucleus are satisfied, while the relativistic effects are taken into account by evaluating the 
matrix elements of the potential in terms of in-medium Dirac spinors. 

In the latter approximation scheme, the Dirac effects are treated directly for the finite nucleus, whereas the correla- 
tion effects are deduced from nuclear matter via a local density approximation. This treatment of the correlation effects 
can be accomplished by defining effective meson exchange interactions, which are obtained from the Dirac-Brueckner 
self-energy components. Therefore, the coupling constants of such an effective interaction are density dependent, since 
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these self-energy components are density dependent [24|. In this way, a semi-phenomenological relativistic density 
functional can be constructed. This approximation scheme to describe finite nuclei will be further explored in this 
paper. 

These semi-phenomenological relativistic density functionals can be divided into two groups, semi-phenomenological 
density dependent relativistic Hartree (DDRH) theories and semi-phenomenological density dependent relativistic 
Hartree-Fock (DDRHF) theories. However, between the semi-phenomenological DDRH theories and the DBHF ap- 
proach two essential differences exist concerning the structure of the self-energy in nuclear matter [TH [l7j . The first 
difference is the absence of a spatial contribution of the vector self-energy Ey in the DDRH theory, which is present in 
the DBHF approach. This self-energy component originates from Fock exchange contributions which are not present 
in the DDRH theory. The second difference is that the DBHF self-energy terms explicitly depend on the particle's 
momentum, a feature which is also absent in the DDRH theory. This momentum dependence reflects the non-locality 
of the DBHF self-energy terms, which originates as well from the Fock exchange terms as from non-localities in the 
underlying effective NN interaction, the G- matrix. In the work of Ref. [TH, HH, H||, renormalizations were intro- 
duced to compensate for these differences. However, a more fundamental solution at the level of the structure of the 
self-energy is preferable and will be investigated in this study. 

An improvement of the simple mean field studies of the DDRH studies is the density dependent relativistic Hartree- 
Fock (DDRHF) theory with its inclusion of the Fock terms, since the essential differences between the DBHF approach 
and the DDRH theory concerning the structure of the self-energy in nuclear matter discussed above may be cured by 
the inclusion of the Fock exchange terms. Due to the presence of Fock terms in the DDRHF theory, one obtains a 
spatial contribution of the vector self-energy Sy and momentum dependent self-energy components as in the DBHF 
approach. However, first attempts of such DDRHF models based on microscopic approaches were not so successful as 
the DDRH models [13, [H| ■ A possible reason is that only the isoscalar coupling functions are density dependent and 
the other coupling functions, such as of the 7r-meson and the p-meson, remain density independent in these DDRHF 
theories [13, HI]. In addition, the 5-meson is absent in these DDRHF theories. However, in DDRH theories this 
coupling provides a mechanism to account for the differences in the scalar self- energies, i.e. in the corresponding 
effective Dirac masses, for neutrons and protons in isospin asymmetric nuclear matter (29|. Therefore, the DDRHF 
theories can be further improved by including the isovector <5-meson and extent the density dependence to all coupling 
functions, which has not been done so far for DDRHF theories based on microscopic calculations. 

The aim of this study is not to provide an "optimal" DDRHF parameterization, which could then be used in studies 
of finite nuclei. Instead we are investigating various models, which are all fitted to reproduce the DBHF results for the 
relativistic components of the nucleon self-energy at the Fermi surface. We than compare these models in predicting 
other observables for nuclear matter and finite nuclei. Therefore it is the aim to explore the limits and the reliability of 
such DDRHF description to provide a reliable approximation scheme for a Dirac-Brueckner-Hartree-Fock description 
of finite nuclei. 

The plan of this paper is as follows. The relativistic DDRHF theory is discussed in SecHH which also introduces the 
various stages of the meson exchange models considered and the resulting parameterization of the density-dependent 
coupling constants. Results for the structure of infinite matter and finite nuclei are presented and discussed in Sec. Mil 
Finally, we end with a summary and the conclusion in Sec. IIVI 



II. DENSITY DEPENDENT RELATIVISTIC HARTREE-FOCK. 



A Lagrangian density of an interacting many-particle system consisting of nucleons and mesons is the starting 
point of a DDRHF theory. The Lagrangian density of the DDRHF theory presented here includes as well the isoscalar 
mesons a and uj as the isovector mesons S and p. Furthermore, one has the pseudo-vector meson 7r. Therefore, the 
Lagrangian density consists of three parts: the free baryon Lagrangian density Cb, the free meson Lagrangian density 
Cm, and the interaction Lagrangian density Ci n t. 

£■ = £-B + Cm + £i n t, (1) 

which takes on the explicit form 

C B = *(i 7M ^ - M)V, (2) 

C M = \ £ (w^-m?*?)-! £ [iF^F^-mlA^A^), (3) 

Ant = -9*^*^ - .g 5 *T* 5 * - -^-*T7b7 M [^* w ]* - g^^A^ - g^ry^A^ (4) 
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with the field strength tensor -F^w = 9^Af K \ v — <9„A( K ) M for the vector mesons. In the above Lagrangian density 
the nucleon field is denoted by ^ and the nucleon rest mass by M. The scalar meson fields are Q a and and the 
vector meson fields are A/ u -s and At p y Furthermore, the pseudo- vector tt field is denoted by <!>„•. Moreover, the bold 
symbols denote vectors in the isospin space acting between the two species of nucleons. The mesons have rest masses 
m a) m u , mg, m p , and and couple to the nucleons with the strength of the coupling constants g a , gs, g u , gp, and 



A. Nuclear Matter 

To obtain the field equations for the nucleons, we minimize the action for variations of the fields $ included in the 
Lagrangian density of Eq. (TTJ 

h r 

dt / d 3 x£(^(x),d^(x),t) =0. (5) 
•> 

Finally the following field equation is obtained 

d dC dC dC5p_ 



= 0. (6) 



Without density dependent meson-baryon vertices, Eq. (j6|) reduces to the normal Euler-Lagrange field equation, 
since the third term vanishes. Evaluating Eq. ([6]) for the field \& , we obtain the Dirac equation. This Dirac equation 
for the nucleon field can be written as 

(i7 M d" - M - (E + E (r) 7o)) * = 0, (7) 

where E is the nucleon self-energy and E^ is the so-called rearrangement contribution to the nucleon self-energy. 
The self-energy E is obtained as well for density dependent as density independent meson-baryon vertices and reads 

E = (g a <P a + gsr&s + ^rjsj^d^} + g^l^A^ 

V \ (8) 

+ 9p-ri^ p} ) ■ 

for the Lagrangian density in Eqs. (HJ-Q. Furthermore, this nucleon self-energy can be split into different parts with 
well-defined behavior under Lorcntz transformations. Because of the requirement of translational and rotational in- 
variance, hermiticity, parity conservation, and time reversal invariance, the most general form of the Lorentz structure 
of the self-energy is 

E(fc) = E s (fc)- 70 E o (fc)+ 7 .kE v (fc), (9) 

where E s , E , and E v components are Lorentz scalar functions. Therefore, it is practical to define the following 
effective quantities 

k* = k(l +%leY± v (k)), (10) 
m*(fc) = M + KeE s (fc), (11) 

and 

E*(k) =E(k) + $ei; o (k). (12) 
In the on-shell case, the effective energy can also be written as 

E* (k) 2 = k* 2 + m* (k) 2 . (13) 
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On the level of the Hartree-Fock approximation, the contributions to the self-energy of the Lagrangian density pre- 
sented in Eqs. (J!])-® are 

£ s ,i(fc) = -( — ) (p s ,n + Ps,p) - [ —) ^2 (p s ,i ~ Ps,j) 



(4 

+(2-%) 



m*(q) 
~E*(q) 



5 i3 [glQ tr (k,q)-4gie u (k,q)] 



^- ) m^e.ffc, ?) - 4 3 2 9 p (fc, g) + 3 2 9 A -(fc, g ) 
m. 



(14) 



and 



% [^e ff (fc,g)+2^e a) (fe,g)] 



+ (2-%) 



(47rfc) 2 ^ 7 E*(q) 



^ [2^ r<T (fc,g)+4 5 2r w (fc,g)] 



+ (2 - ^) 
+2g 2 s T s (k,q) 



- 2 



((fc 2 + q 2 )T n (k, g) - fcg 2 9 ff (fc, g)) + 4 3 2 r p (fc, g) 



The first two terms in Y, s ^ and £ ,z correspond to the Hartree contribution with 



1 



1,3 

3tt 2 



g 2 d<7 



m*(q) 
E*{q) 7 



The remaining expressions are due to the Fock contributions, where the abbreviations 

A(k,q)=q 2 + k 2 -(E(k)-E(q)) 2 , 

( A(k,q) +m 2 + 2kq 



&i(k,q) = In 



\A(k,q) +m 2 - 2kq 



and 



w , , ((A(k,q)+m 2 )e t (k,q)) 

ri(fc,g) = - k q 



(15) 



(16) 

(17) 
(18) 

(19) 
(20) 

(21) 



are used. 

The T,( r ) term in Eq. ([7]) will only be present, if meson-baryon vertices are density dependent, and is generated by 
the third term in Eq. ([5]). This rearrangement contribution reads 



E (r) = (^$ a V + ^* T *,« + -L^Vt^ [8*^]* 
V Op Op Op 



(22) 
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meson i 


m [MeV] 


at hi a di 


a 


550 
782.6 


9.28996 -0.8403 0.25791 -0.015436 
10.23446 0.23622 -0.1503 0.0297512 


a 


550 
782.6 


8.86711 -1.1508 0.39224 -0.034995 
8.37380 0.71233 -0.2986 0.0469882 



TABLE I: Parameter set for the cru) model from the DBHF approach in Ref. [T^]. The upper part of the table contains the 
parameters to be used in the Hartree approach, whereas the lower part refers to the DDRHF model. 



meson i 


m [MeV] 


ai 


hi 


Ci 


di 


a 


550 


8.65683 


-1.0265 


0.35362 


-0.0408585 


LJ 


782.6 


8.62502 


0.53724 


-0.2312 


0.0386465 


TV 


139 


1.00265 









TABLE II: Parameter set for the auiv DDRHF model from the DBHF approach in Ref. [TjJ. 



These rearrangement contributions are essential to provide a symmetry conserving approach, which implies that 
energy-momentum conservation and thermodynamic consistency like the Hugenholtz - van Hove theorem arc 
satisfied [30|. However, these rearrangement terms do not contribute to the energy per nucleon, 



E/A = T + V-M, 



where the kinetic energy per nucleon is 

T = 

and the potential energy per nucleon is 

1 



Mp s 



r 2 

i—n : p " 



k p i 



E*(q) 



q 3 dq 



1 

PB 



V = 



2n 2 p B ^ 

i—n,p 



E*(q) 



T, s ,i(q) - £ Q ,i(g) + 



E*(q) 



qT, Vti (q) q dq 



(23) 



(24) 



(25) 



with p s = p Stn + p s , p and p B = p n + P P - 

Furthermore, the rearrangement terms can be added to the time-like vector self-energy in Eq. (|15[) . Therefore, 
it should be noted that these rearrangement terms are sometimes included in the definition of the time-like vector 
self-energy T, (k). 



B. Parameterization 



We have constructed three different DDRHF models: a <toj model, a auiir model, and a auJirpS model. The effective 
coupling constants for the mesons are determined by requesting that the HF expression for the scalar self-energy S s (/c) 
and the time-like vector self-energy Eo(fc), i.e. without the rearrangement terms E' r \ calculated at the Fermi surface 
reproduce the corresponding results of a DBHF calculation using Bonn A [15j. The reason that the rearrangement 
terms are not included in the fit is that the DBHF approach has no rearrangement contributions. The DBHF data 
at densities of ps = 0.100, 0.197, 0.313, and 0.467 fm -3 are used to obtain the parameters of the various models. The 
density dependent couplings of the auj and the guj-k model are obtained from the DBHF results in isospin symmetric 
nuclear matter (Y p = 0.5), whereas couplings of the aunrpS model are determined from results in isospin asymmetric 
nuclear matter with a proton fraction of Y p = 0.4. For completeness two different Hartree (DDRH) models have been 
constructed in a similar way. 

In order to make these parameterizations easily accessible, we have parameterized the density dependence of the 
coupling constants by 



Sz(rb) = a% + + CiX 2 + diX 3 , 



(26) 
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meson i 


m [MeV] 


a,i 


bi 


a 


di 


a 


550 


9.26408 


-0.79477 


0.24135 


-0.013504 




782.6 


10.18505 


0.31476 


-0.19474 


0.037339 


S 


983 


9.58644 


-4.34658 


1.56009 


-0.208196 


P 


769 


9.51803 


-3.36524 


1.36753 


0.158856 


a 


550 


4.99231 


3.27705 


-1.4929 


0.230977 


U] 


782.6 


0.67868 


10.3068 


-4.6471 


0.701354 


TV 


139 


1.00265 








s 


983 


12.3664 


-8.9781 


3.25210 


-0.426997 


p 


769 


15.2464 


-11.1504 


4.27083 


-0.543179 



TABLE III: Parameter set for the atoirpS model from the DBHF approach in Ref. [la) . The upper part of the table contains 
the parameters to be used in the Hartree approach, whereas the lower part refers to the DDRHF model. Note that the n does 
not contribute in the Hartree approximation. 
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FIG. 1: (Color online) Density dependence of the isoscalar coupling functions of the aco (dotted), aun (dashed), and the 
<7Ci;7rp<5(dashed-dotted) DDRHF model. Left: the a coupling function. Right: the u coupling function. 



with x = pb/po and po — 0.16 fm -3 . The parameters of the coupling functions are fitted, except the one of the tt 
coupling function. The 7r coupling constant is fixed to the free value. In addition, the masses of the mesons are chosen 
to be identical to those of the Bonn A potential. All parameters are summarized in table [I] for the era; model, in table 
[IT] for the aunr model, and in table UTT1 for the auirpS model. 

The density dependence of the isoscalar coupling functions is displayed in Fig. [TJ In the ulj model, the a coupling 
function is decreasing with increasing density, whereas the u coupling function is slightly increasing with increasing 
density. These features can be interpreted in the following way: 

A significant part of the medium-range attraction contained in the Brueckner G-matrix is due to the iterated ir- 
exchange term. Pauli blocking and dispersive corrections of the nucleon-nucleon propagator in the nuclear medium 
yield a quenching of these iterated n— exchange terms with increasing density (see e.g. [l[ and @)- In the auj model 
this medium range attraction is described in terms of the cr-exchange. Therefore, the quenching of the iterated 
7r-exchange terms leads to a reduction of the coupling constant for the a with increasing density. 

The same Pauli blocking effects and dispersive corrections in the NN propagator also reduce the correlation effects 
in the relative wave functions of two nucleons at short distances. Therefore, these short-range correlations are less 
efficient at higher densities to minimize the repulsive short-range components of the bare NN interaction. This leads 
to a slight increase of the effective coupling constant for the to with increasing density. 

Next, we will compare the isoscalar coupling functions of the au> model with that of the aujn model. One finds that 
the inclusion of the ir exchange has a small effect on the isoscalar coupling functions, i.e. the a coupling function is 
slightly decreased, whereas the u coupling function is slightly enhanced compared to the au> model. 

However, in the auiirpS model the density dependence of the isoscalar coupling functions dramatically change 



7 



10 




15 



10 



ocoitpS DDRHF 




0.1 



0.2 

P R [fm 



0.3 



0.4 



FIG. 2: (Color online) Density dependence of the isovector coupling functions of the aunpS DDRHF model. Left: the p 
coupling function. Right: the 8 coupling function. 



compared to the au: model and the aum model. The effective coupling constants for the isovector mesons are rather 
large in particular at small densities as displayed in Fig. [2] This reflects the significant differences between the isospin 
T = 1 and T = interaction, which is required to describe nucleon self-energies in isospin asymmetric nuclear matter. 
The resulting interplay between the exchange terms for isoscalar and isovector mesons spoils the simple picture to 
explain the density dependence of the coupling constants for the scalar mesons discussed above. 



III. RESULTS AND DISCUSSION 



A. Nuclear Matter 



In this section we will present the nuclear matter results of the previously constructed DDRHF models and compare 
them to the original DBHF approach of Ref. [15j . which the various parameterizations attempt to reproduce. First, the 
momentum dependence of the DDRHF self-energy components will be considered. Note that the parameterizations 
are fitted to the self-energy components at the Fermi momentum kp. The momentum dependence of the DDRHF 
self-energies exclusively originates from the various Fock exchange terms, whereas the momentum dependence of the 
original DBHF self-energies is due to Fock exchange terms but also due to the non-locality and energy-dependence of 
the underlying G-matrix. 

This momentum dependence of the self-energy components in isospin symmetric nuclear matter of a density pg = 
0.197 fm -3 for the various DDRHF models is displayed in Fig. [3j The Hartree approximation ignores all Fock exchange 
terms and therefore provides a scalar self-energy which is independent on the nucleon momentum k. Also for the auj 
DDRHF model, the scalar self-energy E s turns out to be almost momentum independent. The inclusion of the tt 
improves the momentum dependence of the scalar self-energy S s . 

However, for the time-like vector self-energy Eo the opposite behavior can be observed, i.e. the inclusion of the tt 
deteriorates the results. Therefore, a awir DDRHF model will in principle not be able to reproduce the momentum 
dependence of the DBHF self-energy components. Although the momentum dependence in the aixnrpd model is in 
better agreement with the DBHF approach than in the other two models, its momentum dependence is only about 
one half of the momentum dependence observed in the DBHF approach. 

As the models are fitted to reproduce the values of the the self-energy components at the Fermi momentum 
k = kp = 1.43 fm~ 3 all models converge at this point. The auJirpS DDRHF model is fitted in isospin asymmetric 
nuclear matter, therefore a small deviation for this DDRHF model can be observed in isospin symmetric nuclear 
matter. 

The other essential difference between the relativistic Hartree model DDRH and the DBHF approach is the absence 
of a spatial contribution of the vector self-energy Ey in the DDRH theory. Such a contribution is present in the 
DDRHF theory due to the Fock exchange terms. In our procedure to determine the various DDRHF models this 
self-energy component is not fitted. Results of various parameterizations are displayed in Fig. S) One observes that 
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k [fm"'] k [fm" 1 ] 

FIG. 3: (Color online) Momentum dependence of the scalar self-energy E s and the time-like vector self-energy Eo for the 
acu (dotted), <jaj7r(dashed), and the cruj7rp<5(dashed-dotted) DDRHF model is plotted at a density of pb = 0.197 fm -3 . In 
addition, the Hartree model with a and ui mesons (dashed-dotted dotted) is given. Furthermore, the corresponding self-energy 
components of the DBHF approach in Ref. t lS] (solid line) are presented in this figure. Left: the scalar self-energy E s . Right: 
the time- like vector self-energy Eo. 
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FIG. 4: (Color online) Density dependence of the spatial vector self-energy E„ for the aoj (dotted), aun (dashed), and the 
cru;7rp<5(dashed-dotted) DDRHF model is plotted in isospin symmetric nuclear matter. In addition, the corresponding self- 
energy component of the DBHF approach in Ref. [1EJ is presented in this figure. 



the models with ir exchange yield a stronger spatial component and show a much better agreement with the DBHF 
result than the aco model. Therefore, it can be concluded that the inclusion of the ir exchange in DDRHF models is 
essential in reproducing the spatial vector self-energy. 

The self-energy components are needed to calculate the energy per nucleon as can be seen from the Eqs. (I25|) - (l?51) . 
The energy per nucleon for isospin symmetric nuclear matter and pure neutron matter are displayed in Fig. [5] The 
various approximation schemes reproduce the DBHF results for symmetric matter reasonably well at medium and 
larger densities. At densities below 0.8 fm -3 all parameterizations tend to underestimate the binding energy per 
nucleon calculated in the DBHF approach. The au> model is less attractive than the DBHF and the other approaches 
at larger densities. This is partly due to the weak spatial vector component but also due to the fact that the difference 
between the self-energy components E s and Sq tends to be too small at low momenta in this approach as can be seen 
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p B [fm" 3 ] p B [fm' 3 ] 

FIG. 5: (Color online) Energy per Nucleon is plotted as a function of the density for the cjuj (dotted), ausn (dashed), and the 
aa;7rp5(dashed-dotted) DDRHF model. In addition, the energy per nucleon of the DBHF approach in Ref. [15| (solid line) is 
given. Left: isospin symmetric nuclear matter. Right: pure neutron matter. 





Hartree 
auj aunvpS 


DDRHF 


DBHF 


po [fm" 3 ] 
E/A [MeV] 


0.1572 0.1613 
-14.44 -14.76 


0.1508 0.1623 0.1583 
-13.56 -15.25 -16.34 


0.181 
-16.15 



TABLE IV: Saturation properties of nuclear matter for the various models compared to the DBHF results of Ref. [l4, 15]. The 
quantities listed include the saturation density po and the binding energy E/A at saturation density. 

in Fig. O 

Therefore, at saturation density nuclear matter is too weakly bound as can be seen in Table. IIVI In addition, the 
saturation density is shifted to lower densities compared to the original DBHF results. This shift of the saturation 
density to a lower density than in the DBHF approach is also observed in the other DDRHF models. 

For pure neutron matter, one finds in Fig. [5] that the isovector mesons p and S are important to describe the DBHF 
results at low densities in particular. The influence of these mesons gets smaller at high densities, which is already 
indicated by the density-dependence of the coupling constants shown in Fig. [2] 

B. Finite Nuclei 

For the study of finite nuclei, we account for the density dependent correlation effects in a relativistic HF calculation 
by employing the coupling constants calculated at the local density. The density profile ps(r) is determined from the 
result of the relativistic HF calculation in a self-consistent manner (27j . Furthermore, the rearrangement self-energy 
contribution is taken into account, since it is important to get appropriate single-particle energies and wave functions 
in finite nuclei. To solve the Dirac equation for finite nuclei in coordinate space, the radial functions are expanded 
in a discrete basis of spherical Bessel functions. This discrete basis forms a complete orthonormal basis in a sphere 
of radius D. This radius D is chosen to be 30 fm, which is large enough to ensure that the results for the bound 
single-particle states are independent on D. 

We have investigated some closed-shell nuclei with our DDRHF model. The results obtained for the binding energy 
and charge radius are presented in Table |Vj The center of mass correction is included in the displayed binding energies 
and the charge radii have been evaluated from the proton density, assuming a radius of 0.8 fm for the charge radius 
of the proton. 

The au DDRHF model yields with too little binding for all nuclei considered. This outcome is expected from the 
nuclear matter results, since this model is too repulsive in symmetric nuclear matter as can be seen in Fig. [5] The 
explicit inclusion of the 7r-exchange improves the results substantially. The aujirpS model yields too much binding 
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Nucleus 


Hartree 


Hartree 


DDRHF 


DDRHF 


DDRHF 


Exp. 




<JLO 


aunpS 




auin 


aujirpS 






E/A r C n 


E/A r CD 


E/A r C D 


E/A r C D 


E/A r C D 


E/A r C D 




-6.49 2.75 


-6.65 2.72 


-6.00 2.77 


-7.63 2.67 


-8.91 2.62 


-7.98 2.74 


40 Ca 


-7.10 3.47 


-7.29 3.43 


-6.54 3.49 


-8.08 3.39 


-9.18 3.34 


-8.55 3.48 


4S Ca 


-7.29 3.51 


-7.20 3.47 


-6.64 3.53 


-7.67 3.43 


-8.81 3.40 


-8.67 3.47 


90 Zr 


-7.20 4.29 


-7.26 4.23 


-6.56 4.32 


-7.70 4.20 


-8.83 4.16 


-8.71 4.27 



TABLE V: Results for the binding energies per nucleon in MeV and the charge radii in fm for the various DDRH and DDRHF 
models. The calculated energies have been corrected by subtracting the spurious energy of the center of mass motion. The 
experimental values are taken from Ref. [25l ]. 





(JU1 




auin 




auinpS 




exp. 




Orbital 


Proton 


Neutron 


Proton 


Neutron 


Proton 


Neutron 


Proton 


Neutron 


160 










Osi ii 


-27.62 


-31.86 


-34.64 


-39.05 


-37.06 


-41.50 


-44±7 


-47 




-13.16 


-17.09 


-16.55 


-20.65 


-18.73 


-22.90 


-18.451 


-21.839 


0pi/2 


-9.35 


-13.20 


-12.75 


-16.78 


-14.23 


-18.34 


-12.127 


-15.663 


40 Ca 










0Si/2 


-33.13 


-41.45 


-41.28 


-49.91 


-43.15 


-51.77 


-49.1±12 




0P3/2 


-22.11 


-30.06 


-27.79 


-36.00 


-30.19 


-38.43 


-33.3±6.5 




Opi/2 


-19.55 


-27.45 


-25.30 


-33.47 


-27.42 


-35.61 


-32±4 




0d 5 /2 


-10.71 


-18.29 


-13.62 


-21.45 


-15.72 


-23.64 


-14.2±2.5 


-21.30 


lSl/2 


-6.70 


-14.23 


-8.95 


-16.78 


-9.31 


-17.33 


-10.850 


-18.10 


0^3/2 


-6.61 


-14.08 


-9.50 


-17.23 


-11.01 


-18.85 


-8.325 


-15.641 


4S Ca 










Osi/2 


-37.17 


-42.35 


-44.52 


-50.90 


-48.38 


-52.64 






0P3/2 


-27.25 


-31.13 


-32.05 


-36.84 


-36.75 


-39.03 






Opi/2 


-25.18 


-29.51 


-29.80 


-36.46 


-34.10 


-38.50 






0d 5 /2 


-16.32 


-19.48 


-18.38 


-22.09 


-23.22 


-24.10 


-20±1 




lSl/2 


-11.26 


-15.57 


-13.03 


-18.45 


-15.84 


-19.08 


-15.8 




0d s/2 


-12.55 


-16.18 


-14.47 


-20.93 


-18.68 


-22.44 


-15.3 




0/7/2 


(-4.90) 


-7.89 


(-4.55) 


-7.82 


(-8.75) 


-8.95 







TABLE VI: Single particle energies for the orbital levels of O, Ca, and Ca nucleus derived from RHF calculations using 
various models. The experimental values are taken from Ref. [27l. l3lj. 



energy for the light nuclei, whereas a good agreement to the empirical values is obtained for the heavier nuclei. 
Furthermore, all models predict charge radii, which are in good agreement with experimental values. 

Single particle energies for the orbital levels of the le O, 40 Ca, and 48 Ca nucleus are presented in Table IVII for 
the various models. The aui DDRHF model predicts absolute values of single-particle energy, which are too small. 
The single particle energies of the auin and auinpS model are in good agreement with experimental data. They are 
in better agreement with experimental values than the energy levels calculated in Refs. [27l . l28j . which lie to deep. 
However, the order of the ls-shell and Od-shell deviates from experiment. The spin-orbit splitting plays a key role in 
the ordering of these shells. The spin-orbit splittings of some nuclei are given in Table IVIII It can be seen that in 
all the presented DDRHF models the spin-orbit splittings in le O and 40 Ca are smaller than the experimental values. 
However, the agreement is better in 4S Ca and in 90 Zr, in particular for the auinpS model. In works of Ref. [27t l28j. 
it was reported that the n exchange reduced the spin-orbit splitting. We find as can be seen in Table IVIII that the 
inclusion of the n has no noticeable effect on the spin-orbit splitting in light symmetric nuclei, whereas in heavier 
asymmetric nuclei a small increase in the spin-orbit splitting can be observed. However, the inclusion of the isovector 
p and S meson leads to a clear increase of the spin-orbit splitting in as well light as heavier nuclei. 

Fig. [5] shows the proton density distribution of the 16 O nucleus. It is found that in the aui model the proton density 
in the interior is smaller than in the other two models and has a longer tail. Therefore, it has a larger charge radius as 
can be seen in Table IVl Furthermore, the smaller maximum density in the aui model compared to the other models 
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Nucleus 


aw 


(TUJTV 


aujnpS 


exp. 




3.81 


3.8 


4.5 


6.3 


40 Ca 


4.1 


4.12 


4.71 


7.2 


4S Ca 


3.77 


3.91 


4.54 


4.3 


90 Zr 


1.03 


1.22 


1.51 


1.5 



TABLE VII: Spin-orbit splittings of protons for the Op-shell in u '0, the Od-shell in 40 Ca and 48 C"a, and the lp-shell in 90 Zr. 
The experimental values are taken from Ref. [28l ]. 
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FIG. 6: (Color online) Proton density distribution of ie O nucleus for the auj (dotted), aunv (dashed), and the CTCJ7rp5(dashed- 
dotted) DDRHF model. 



can be expected from its smaller saturation density in nuclear matter. 



IV. SUMMARY AND CONCLUSION 



Although Dirac-Brueckner-Hartree-Fock (DBHF) calculations have been quite successful in describing nuclear mat- 
ter, corresponding DBHF calculations have not yet been performed for finite nuclei as such calculations still seem to 
be too complex. Various attempts have already been made to approximate such DBHF calculations by employing 
some kind of local density approximation by parameterizing the results of DBHF calculations of nuclear matter in 
terms of either a Dirac Hartree (DDRH) model or a Dirac Hartree Fock (DDRHF) model with density-dependent 
coupling constants for various mesons considered. 

It has been the aim of this investigation to compare these different models. The parameters of these DDRH and 
DDRHF models were all fixed to reproduce at each density the same scalar E s and time-like vector component Eo 
component of the self-energy obtained in the DBHF calculation for nucleons with momentum equal to the Fermi 
momentum kp. 

While the DDRH approach yields no momentum dependence of the self-energies, the DDRHF model reproduces the 
qualitative features of the momentum dependence obtained in the DBHF calculations. However, only 50 percent of this 
momentum dependence can be related to the Fock exchange terms. The remaining part of the momentum dependence 
in DBHF is due to the non-locality of the DBHF G-matrix. It remains as a challenge for further investigations to 
account for this non-locality in such a way that the complete energy- and momentum-dependence of the self-energy 
is described in nuclear matter using a representation, which can be transferred to calculations of finite nuclei. 

The spatial vector component of the nucleon self-energy E w of the DBHF is dominated by the 7r-exchange term and 
therefore can be accounted for in DDRHF models, which include the it explicitly. The exchange of isovector mesons 
p and 5 is important to describe the differences of correlation effects in isospin T — and T = 1 channels. Therefore 
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a DDRHF model accounting for a, co, tt, p and S exchange is required to reproduce the features of DBHF calculations 
for symmetric and asymmetric nuclear matter. 

Fair agreement with empirical data is obtained when this parameterization of the full aunrpS model is employed in 
Dirac Hartree Fock calculations of finite nuclei with density dependent coupling constants for light and heavy nuclei. 
The parameterization of the coupling constants is presented in form which makes it accessible for other users. 
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